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Abstract In California, severe precipitation events 
(SPEs) are often associated with winter season 
atmospheric rivers. These SPEs can generate hurri- 
cane-scale precipitation, creating a variety of natural 
hazards such as floods and landslides. For California, 
the most complete SPE-flood record yet generated is 
a 9000-year paleoflood reconstruction from the Santa 
Barbara Basin (Du et al. in Mar Geol 397:29-42, 
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2018). Finding terrestrial counterparts to compare 
to the Santa Barbara Basin is a challenge in South- 
ern California where lake basins are rare and the 
terrestrial spatiotemporal signature of SPE floods 
remains largely unconstrained. Here, we present five 
1—2 m sediment cores from a sag pond along the San 
Andreas Fault in the southern Carrizo Plain. The 
most complete record (core CLPC21-4) was selected 
for a variety of chronological, sedimentological, and 
biological analyses. Principal Component Analy- 
sis on these data reveals two end member sediment 
unit types: event versus ambient sedimentation. We 
focus on the two thickest sediment units likely gener- 
ated by SPEs. These two units (EU 1 and EU 2) are 
especially distinct within CLPC21-4, characterized 
by erosive bases, above average sand content, normal 
grading, low magnetic susceptibility, and low total 
organic matter. Moreover, they are visually apparent 
across all five cores within the sag pond and thin from 
source to sink. Age control for the two units is con- 
strained by AMS |4C dates on discrete organic mate- 
rials and supported by the identification of Erodium’s 
first appearance ca. 1750-1765 CE in the Santa Bar- 
bara region. Using these age constraints, we infer the 
maximum limiting age range for Event Unit 1 from 
1470 to 1640 CE and Event Unit 2 from 1740 to 1800 
CE. Within the limits of dating, site-specific proxy 
sensitivities, and inherent meteorological heterogene- 
ity, we propose a probable correlation to Santa Bar- 
bara Basin flood events at 1525 CE and 1760 CE. Our 
results suggest that sag ponds may represent a viable 
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and untapped paleoclimatic archive for California. 
Future work will focus on a latitudinal series of sag 
ponds to determine the spatiotemporal sequence and 
correlativity of SPEs in the sediment record. 


Keywords Grain size - Sedimentology - Event 
deposition 


Introduction and background 


When you think of natural hazards and California, 
earthquakes dominate the narrative. But in terms of 
socioeconomic impact, floods are the single greatest 
natural hazard threat to California (Dettinger 2011; 
Porter et al. 2011; Corringham et al. 2022; Huang 
and Swain 2022). For example, in 1861-1862 AD, 
California experienced an historically unprecedented 
series of winter storms generating extreme pre- 
cipitation and statewide flooding (Engstrom 1996). 
This 45-day period of record storms flooded nearly 
every major California metropolitan area at the time. 
Moreover, now dry lakes in the Mojave re-emerged, 
Tulare Lake expanded, and significant coastal ero- 
sion occurred along the California Coast (Porter 
et al. 2011; Reynolds et al. 2018). The USGS mod- 
eled this scenario (i.e., the ARkStorm, or ~ 1000-year 
event) as if it happened today and estimated the pro- 
jected cost approximately 2xthe estimated impact 
of the potential “big” San Andreas Fault earthquake 
(Porter et al. 2011). Future climate projections sug- 
gest that extreme precipitation accumulation amounts 
may increase by 10-40%, extreme precipitation vari- 
ability will increase, and extreme events like that in 
1861-1862 AD may increase in frequency by up to 
three-fold during the twenty-first century (Polade 
et al. 2017; Swain et al. 2018; Gershunov et al. 2019; 
Huang et al. 2020; Huang and Swain 2022). 

At present, our knowledge of California’s flood 
history is derived from marine or coastal sediments 
(Fig. 1). Several important records illustrate flood 
events such as a 9000-year flood event history from 
Santa Barbara Basin (SBB; Schimmelmann et al. 
1998; Hendy et al. 2013, 2015; Du et al. 2018; Sarno 
et al. 2020), a 7000-year potentially flood-related tur- 
bidite sequence in the Santa Monica Basin (Romans 
et al. 2009), a mid-Holocene flood event layer in Lake 
Elsinore (Kirby et al. 2021), and evidence for the 
1861-1862 AD over wash deposit from Carpinteria 
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Marsh (south of Santa Barbara, CA; Reynolds et al. 
2018). Du et al. (2018) show that flood frequency 
and magnitude are in general positively correlated 
throughout the Holocene. Moreover, the frequency 
of floods tends to cluster and shows regional corre- 
lations to pluvials from terrestrial archives (Kirby 
et al. 2010, 2012, 2014, 2015; Hiner et al. 2016). 
More recently, Du et al. (2021) revealed that the SBB 
flood record shows a strong Holocene-scale connec- 
tion to El Nifio-Southern Oscillation (ENSO) vari- 
ability and Pacific ocean—atmosphere dynamics. This 
relationship agrees with lower-resolution terrestrial 
records from California that show Holocene-scale 
couplings between ENSO and Pacific ocean—atmos- 
phere dynamics for droughts and pluvials (Masters 
2006; Covault et al. 2010; Miller et al. 2010; Kirby 
et al. 2014, 2015, 2023; Hiner et al. 2016; MacDon- 
ald et al. 2016). 

Perhaps most revealing from the SBB flood record 
is the relatively subtle sediment signature associated 
with the 1861-1862 AD “ARkStorm-type” event. If 
event thickness is an approximate indicator of event 
magnitude (reasonable as a first approximation; Saba- 
tier et al. 2022), the ARkStorm scenario is likely 
underestimating the scale of future SPEs and their 
associated flood risk in California. In fact, the SBB 
flood record is characterized by at least 30 flood- 
event layers that are thicker than the rather undersized 
1861-1862 AD event (Du et al. 2018). This realiza- 
tion represents a potentially important mismatch 
between the geologically observed record and mod- 
eling output. Critical to predicting future flood sce- 
narios is the ground-truthing of predictions and mod- 
els in the geological record. We propose to address 
this very issue using sediments from a sag pond along 
the San Andreas Fault in the southern Carrizo Plain 
National Monument (Fig. 1). Moreover, despite a 
very well dated marine record offshore from Santa 
Barbara, the terrestrial spatiotemporal signature of 
SPE floods remains largely unconstrained. 


Regional climatology 


The climate of Southern California is Mediter- 
ranean—cool, wet winters and hot, dry summers 
(Bailey 1966). Year-to-year winter precipitation 
variability is related to tropical and north Pacific 
ocean—atmosphere dynamics and their modulation of 
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Fig. 1 a Regional and b local maps with core site locations 
(Google Earth Pro v.7.3.6.9345, 4/22/21 image date) and c 
major drainage area noted (see Supplementary material 1). a 
CPSP =Carrizo Plains Sag Pond, ZL=Zaca Lake (Kirby et al. 
2014), AL= Abbott Lake (Hiner et al. 2016), SBB = Santa Bar- 


the winter-storm track over the western US (Cayan 
and Peterson 1989; Castello and Shelton 2004; Han- 
son et al. 2006). For example, both ENSO and the 
Pacific Decadal Oscillation modulate the mean win- 
ter position of the eastern Pacific subtropical high and 
thus the average latitude of the Pacific winter storm 
track (Cayan and Peterson 1989; Castello and Shelton 
2004; Hanson et al. 2006; Wise 2010, 2016). When 
conditions force the storm tracks further south, they 
produce higher-than-average precipitation and greater 
river discharge in Southern California, whereas a 
more northerly position results in lower-than-average 
precipitation and less river discharge (Cayan et al. 
1998; Hanson et al. 2006; Gray et al. 2015). 

In California, nearly all major floods in the twen- 
tieth century are attributed to winter season atmos- 
pheric river (AR) storms (Dettinger 2011; Dettinger 
et al. 2011; Ralph and Dettinger 2011; Konrad and 
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bara Basin (Sarno et al. 2020), SG=San Gabriel regional tree 
rings (MacDonald et al. 2008; Meko et al. 2017), LE=Lake 
Elsinore (Kirby et al. 2010, 2018), ML= Mojave Lakes (Enzel 
et al. 1992; Miller et al. 2010; Honke et al. 2019) 


Dettinger 2017; Corringham et al. 2019). These nar- 
row, moisture-rich plumes transport enormous quan- 
tities of water vapor from the tropics and subtropics 
into the mid-latitudes (Ralph and Dettinger 2011). 
ARs converge with coastal ranges, experience oro- 
graphic lifting, and generate precipitation amounts 
several standard deviations above normal, often lead- 
ing to landslides, debris flows, floods, and possibly 
incising new stream channels (Grant-Ludwig et al. 
2010). Only land-falling hurricanes of the North 
Atlantic generate comparable precipitation amounts 
in North America (Dettinger et al. 2011). Although 
much less frequent than winter-season ARs, dissipat- 
ing eastern tropical Pacific cyclones, such as Hur- 
ricane Kathleen in 1976, may impact the study area 
(Corbosiero et al. 2009). Currently, it is not possible 
to unequivocally differentiate the precipitation source 
of paleo-floods; however, statistically a winter-season 
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AR is more likely the cause of flooding at the study 
site than is a late-summer tropical cyclone. 


Study site 


The study site is a sag pond (hereafter Carrizo Plain 
sag pond; CPSP: N34°59'30.35"; W119°29'0.13") 
located along the San Andreas Fault (SAF) between 
the Temblor Range and the Caliente Range at 804 
mas! (Fig. 1). Sag ponds are topographic depressions 
formed along strike-slip faults; they are common 
along the SAF in California. Previous earthquake- 
interval research done on the SAF 36 km NW of 
our study site indicates that the recent average inter- 
val of rupture events is 88+41 years with a slip rate 
of ~35 mm year™! (Akgiz et al. 2010; Grant-Ludwig 
et al. 2010). Rupture events can affect the sedimen- 
tology of sag pond sediments and will be considered 
when interpreting the site’s sedimentology. Although 
the current lake surface area is small (~0.021 km’), 
the sag pond receives runoff from a large drainage 
basin (~9.2 km?; (see Supplementary material 1). 
CPSP is also a terminal basin with no outlets. The 
drainage basin is highly asymmetric with most of the 
drainage originating from the south, west, and north- 
west sides of the lake (see Supplementary material 1). 
Deeply incised channels characterize the south, west, 
and northwest drainage basin. A fan delta-like fea- 
ture is apparent along the lake’s northwestern edge, 
connecting to a visually distinct channel despite con- 
struction of a road bisecting the channel (see Sup- 
plementary material 1). Although observational, the 
incised channels and fan delta-like feature suggest 
periods of active flow in recent times. The eastern 
side of the lake is steep and deeply incised, but it rep- 
resents a much smaller drainage area. Google Earth 
historical imagery from 1985 to 2022 reveal instances 
of standing water in the southeastern sector (1.e., 
modern depocenter) of the lake as recently as Febru- 
ary 2020. Notably, there are only two images from 
1985 to 2002, so the breadth of the lake’s modern 
hydrological fluctuations is poorly constrained. The 
local modern vegetation is predominantly a shrub and 
grassland community (Buck-Diaz and Evans 2011). 
Exotic grass species such as Erodium are common. 
There are no climate data from the study site. As 
a result, we infer its climate from historical mete- 
orological data from the Western Regional Climate 
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Center database using Taft, CA, as our comparison 
site (NWS Cooperative Network Number 048752, 
1948-2008 AD, 17 km NNE of CPSP). Summer 
(JJA) temperatures average 26.7 °C; winter (DJF) 
temperatures average 10.4 °C. Average annual pre- 
cipitation averages 137 mm with 105 mm (or 76.8%) 
of the total occurring between November and March. 


Methods 


Five sediment cores—CLPC21-1 to CLPC21-5— 
ranging from 125 to 156 cm in length were taken 
along the long axis of CPSP using a post-driver 
attached to a 4” diameter aluminum barrel (Fig. 1, 
Table 1). The cores were transported back to Califor- 
nia State University, Fullerton where they were split 
lengthwise, visually described, and digitally pho- 
tographed. Of the five sediment cores, we focused 
on CLPC21-04 (156 cm) for detailed analyses and 
age control, being the most complete core and with 
the thickest event units. In two of the five cores 
(CLPC21-01, -02), the core sediment separated dur- 
ing core extrusion (i.e., noted as core breaks in the 
photographs). Because the core is a single drive, we 
know that there is no missing sediment associated 
with the core break as the break occurs mid-core. 
Although the research focused on CLPC21-4, all five 
cores were examined using visual stratigraphy, mag- 
netic susceptibility, and percent water content. 

Age control for this study was determined on core 
CLPC21-4 using radiocarbon dating (accelerator 
mass spectrometry (AMS) '4C dating) on seventeen 
discrete microscopic and macroscopic (> 125 pm) 
organic materials (e.g., charcoal, seeds, insect parts). 
Samples were measured on the insoluble fraction 
for 4C at the University of California, Irvine Keck 


Table 1 Study site core locations 


Core ID Latitude Longitude 

CLPC21-1 34 59'27.3"N 119 28'56.2W 
CLPC21-2 34 59'27.4"N 119 28'56.1W 
CLPC21-3 34 59'28.9"N 119 28'58.3W 
CLPC21-4 34 59'31.0"N 119 29'00.4W 
CLPC21-5 34 59'32.5"N 119 29'02.5W 


Note: cores 1 and 2 were taken adjacent to one another within 
2-m horizontal distance. Carrizo Lake Cores (804 m ams!) 
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Carbon Cycle AMS Facility (Table 2). All samples 
were treated with an acid—base-acid protocol (IN HCl 
and IN NaOH, 75 °C) prior to combustion. Radiocar- 
bon concentrations are given as fractions of the mod- 
ern standard, D'4C, and conventional radiocarbon 
age, following the conventions of Stuiver and Polach 
(1977). Sample preparation backgrounds have been 
subtracted, based on measurements of 4C_free wood. 
All results have been corrected for isotopic fractiona- 
tion according to the conventions of Stuiver and Pol- 
ach (1977), with 8'°C values measured on prepared 
graphite using the AMS spectrometer. Radiocarbon 
ages were converted to age in calendar years before 
present and age in Common Era (CE) using OxCal 
v.4.4.4 (Ramsey 2021) and the IntCal20 Northern 
Hemisphere radiocarbon age calibration curve from 
Reimer et al. (2020; Table 2). As our focus is on CE 
age, we also report the CE median, minimum, and 
maximum with the latter two ages reported at the 
95.4% range (Table 2). We selected the range of CE 
years (min to max) for any given “C age using the 
highest percent value reported under the 95.4% prob- 
ability distribution curve. '°’Cs was also measured 
in the upper 50 cm of core CLPC21-4 (Table 2). 
Finally, Evodium pollen were identified (presence vs 
absence) at various depth—based on initial '4C dates 
to constrain their likely occurrence—at 10 depths 
from non-event (i.e., ambient) sediment units on core 
CLPC21-4 to determine the first appearance of the 
exotic grass for use as an independent age marker and 
to assess the 4C ages (Table 2; Mensing and Byrne 
1998; Kirby et al. 2014). All 1 cm? pollen samples 
(n= 10) underwent standard pollen-extraction proce- 
dures (Faegri et al. 1989; Bennett and Willis 2001). 
Tracer grains of Lycopodium were added during 10% 
HCl digestion (Batch #1031, Lund University). Pol- 
len identifications were based on modern reference 
collections. Erodium grains were identified only to 
genus. Other pollen taxa present in the samples were 
noted, but not counted. Pollen preservation was poor 
in all samples, though only one sample was unusable. 

Magnetic susceptibility (MS; x 10-7 m? kg~!) was 
measured at | cm contiguous intervals on sediment- 
filled individual 8 cm? plastic cubes using a Barting- 
ton MS2 magnetic susceptibility meter on all five 
sediment cores. The same cubed samples were used 
to determine percent water content and percent dry 
bulk density. Using core CLPC21-4, 1-cm contigu- 
ous samples were combusted at 550 °C and 950 °C 


for two hours to calculate the percent total organic 
matter (% TOM) and percent total carbonate (% TC), 
respectively (Dean 1974). Charcoal per 1-g dry sedi- 
ment were counted on CLPC21-4 from 30 to 156 cm. 
The upper 30 cm were not counted for charcoal due to 
modern, near surface mud cracks that likely compro- 
mise the accuracy of the more-recent charcoal counts. 

Grain size was measured at 1-cm contiguous inter- 
vals on CLPC21-4 following standard pretreatment 
protocols: 30-50 mL of 30% H,O0,, 10 mL of IN HCl, 
and 10 mL of IN NaOH (Leidelmeijer et al. 2021). 
Grain size was determined using a Malvern Master- 
sizer 2000 grain size analyzer attached to a Hydro 
2000G dispersion unit. A ten second sonication pro- 
ceeded each analysis in the Hydro 2000G dispersion 
unit prior to analysis. A silica carbide polishing pow- 
der standard was run twice at the beginning of each 
day, once every ten samples, and once at the end of 
every day to evaluate the equipment’s analytical sta- 
bility and reproducibility over time (n=4127, aver- 
age= 13.11 pm, standard deviation=0.10 pm). All 
data are reported as volume percent and divided into 
10 grain size intervals according to the Wentworth 
scale (Wentworth 1922) as well as d0.5 (0.5 = mean). 

We utilized Principal Component Analysis (PCA) 
based on normalized data and Euclidean distance to 
explore the relationship between magnetic suscepti- 
bility, percent total organic matter, charcoal counts, 
and the ten grain size divisions (n=124 depths of 
shared data) on core CLPC21-4. To help confirm the 
statistical significance of the relationships between 
samples, we calculated a similarity profile (SIM- 
PROF) permutation test. SIMPROF analysis ena- 
bles us to test for structure in multivariate data and 
returns a p-value to determine whether the multivari- 
ate structure manifest in a group of samples are more 
or less similar to each other than would be expected 
if the data were random and lacked structure. A 
more detailed review of this technique can be found 
in Clarke et al. (2008) and Somerfield and Clarke 
(2013). All statistical analyses were conducted using 
PRIMER V7 (PRIMER-E, Plymouth, UK). 


Results 
Visually, all five cores are characterized by relatively 


massive, homogenous sediment (Fig. 2). Interrupting 
this homogeneity are several units that are visually 
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(e.g., darker), sedimentologically (e.g., coarser), and 
structurally (e.g., thicker and with sharp lower bound- 
aries) distinct. Most notable are two units near the 
bottom of each core—hereafter, Event Unit 1 (EU 1: 
151—141.5 cm) and Event Unit 2 (EU 2: 139-128 cm; 
Fig. 2). From north to south in the basin, EU | and 
EU 2 thin and become less coarse. There are other 
visually distinct units nearer the mid-to-top of each 
core, although much thinner than EU 1 and EU 2 and 
not as easy to correlate across cores. Finally, each 
core shows a darker, massive unit of varying thick- 
ness above EU 2 (Fig. 2). To avoid subjective inter- 
pretation of the PCA data, we focus on the two bot- 
tom units only (i.e., the most extreme negative PC1 
values (EU 1 and EU 2)), which are visually corre- 
lated across the core sites. Moreover, we limit our 
detailed sedimentological and chronological efforts 
on the longest core (CLPC21-4), which is the only 
core to capture EU 1 and EU 2 in their entirety and 
with thick enough bedding to extract meaningful 
sedimentological information (i.e.,>1 cm or>one 
point of data; Fig. 3). Cores CLPC21-1, -2 also cap- 
ture both EU 1 and EU 2 in their entirety; however, 
they are much thinner and cannot provide multi-depth 
sedimentological analyses. 

Figure 3 shows all the sediment data including 
the pertinent age control data for core CLPC21-4. 
To evaluate these data objectively and statistically 
in terms of sedimentological characteristics, we uti- 
lized Principal Component Analysis (PCA) based on 
normalized data and Euclidean distance to explore 
the relationships between magnetic susceptibility, 
percent total organic matter, charcoal counts, and 
the 10-grain size divisions (n= 124 depths of shared 
data; see Supplementary material 2). The top 30 cm 
of data was not included in the PCA as it did not 
include charcoal data. PC1 axis accounts for 52.7% 
of the variation within the samples. The PCA plot 
shows that the positive PC1 axis is associated with 
MS (p=1.8, o=0.3), % TOM (p=8.80, o= 1.49), % 
clay (u¥=17.13, o=3.26), % very fine silt (u=24.03, 
o=3.15), % fine silt (u=24.90, o=3.54), and % 
medium silt (u= 19.79, o=3.52); whereas, the nega- 
tive PC1 axis is associated with charcoal (u=23.79, 
o=43.56), % coarse silt (u=18.91, o=6.72% 
very fine sand (u=19.78, o=10.59), % fine sand 
(H=10.12, o=8.44), % medium sand (p=4.80, 
o=9.46), % coarse sand (u=2.71, o=7.07), and 
% very coarse sand (u=0.44, o=1.98). PC2 is not 
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addressed as it accounts for only 19.4% of the vari- 
ance, and it did not add significant (or meaningful) 
information to our interpretation. 

Plotting PC1 shows seven units where the PC val- 
ues are less than zero being characterized by coarser 
sediment particularly coarse silt, very fine, and fine 
sand. As the PC1 values become more negative, the 
content of coarse and very coarse sand and the occur- 
rence of charcoal both increase; although, charcoal is 
generally found only at the top of the most PC1-nega- 
tive units. As an additional analysis of the differences 
between the positive and negative PC1 data, we plot- 
ted percent clay versus the five sand classifications as 
well as percent total sand (Fig. 4). For each compari- 
son, there is a notable difference between clay con- 
tent and sand sizes for the PC1 negative versus PC1 
positive value sediments. PC1 negative values are 
consistently coarser in all sand-size fractions than are 
PC1 positive value sediments. Although four negative 
PC1 units are especially prominent in core CLPC21-4 
(i.e., 151-141.5 cm, 139-128 cm, 75-68 cm, and 
48-42 cm), we focus on the two lower units (EU 1 
and EU 2) because they are characterized by the most 
negative PC1 values, they are sedimentologically and 
visually distinct, they are visually correlated across 
all five core site locations, and they are thick enough 
to extract meaningful sedimentological information 
(i.e., >1 cm or>one point of data per event unit). 
This conservative approach to interpreting the PCA 
results reduces interpretative subjectivity and over- 
interpretation of the less prominent PC1-negative 
units. 

Both EU 1 and EU 2 are characterized by 
decreases in magnetic susceptibility, total organic 
matter, and percent clay, and increases in percent 
very fine and fine sand (Fig. 5). However, EU 1 also 
shows an increase in percent medium, coarse, and 
very coarse sand; whereas, EU 2 shows only a small 
increase in percent medium sand and an absence of 
coarse and very coarse sand. Overall, EU 1 is con- 
siderably coarser than EU 2; however, EU 1 also 
contains a spike in percent clay and magnetic sus- 
ceptibility mid-event at 146-145 cm, apparent as a 
visually lighter colored layer in the digital photograph 
(Fig. 5). Structurally, both EU 1 and EU 2 show sharp 
lower boundaries with normal grading to finer sedi- 
ment up-core to the event unit’s upper boundary. EU 
2 also contains well-defined-to-faint laminae between 
136 and 130 cm; similar laminae are not observed 


in EU | except for the thin clay layer noted above. 
Finally, charcoal content increases slightly at the top 
of EU 1, while it rises sharply to the highest values in 
the core in the upper section of EU 2. 

All age control data are shown on Fig. 6 and 
Table 2 and Supplementary material 3. The surface 
age was assigned modern or 2021 CE. The peak in 
'57Cs occurred at 44-43 cm, which was assigned an 
age of 1963 CE. We used the radiocarbon results 
to identify the zone in which Erodium should first 
appear (see Supplementary material 3). Erodium’s 
first appearance is at 108-107 cm and was assigned 
an age of 1755 CE (range: 1750-1765 CE) based 
on Mensing and Byrne (1998) for the Santa Barbara 
region. Finally, an age model was not constructed for 
core CLPC21-4 because our focus is on the timing of 
EU 1 and EU 2 and not the core’s complete deposi- 
tional history. 


Discussion 


Statistics and depositional interpretations for EU | 
and EU 2 


We focus on the grain size, specifically, to interpret 
these depositional environments as changes in grain 
size in lacustrine environments reflect well-estab- 
lished water energy-grain size dynamics such as 
changes in water depth, distance from the shore zone, 
and grain size-energy dynamics associated with run- 
off from the drainage basin (Anderson 1977; Davis 
and Ford 1982; Hakanson and Jansson 1983; Hilton 
1985; Dearing 1991, 1997; Blais and Kalff 1995; 
Kirby et al. 2010, 2015, 2018; Bird et al. 2017; Pribyl 
and Shuman 2017; Shuman and Serravezza 2017). 
In general, finer grain sizes are deposited in lower 
energy environments, such as a lake’s depocenter or 
its most distal point from the shoreline or sediment 
source; whereas, coarser sediments are generally 
associated with higher energy shoreline environments 
or event deposition caused by high energy processes 
such as precipitation-related runoff events (Anderson 
1977; Dearing 1997; Brown et al. 2002; Kirby et al. 
2010, 2014, 2015, 2021, 2023; Pribyl and Shuman 
2017; Sabatier et al. 2022). 

As previously noted, PC1 axis accounts for 52.7% 
of the variation within the samples with the positive 
PC1 axis associated with above average MS, % TOM, 
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<Fig. 2 All five core photographs with Event Units 1 and 2 
highlighted (i.e., thickness of highlight=thickness of event 
unit). Core depth is shown at the bottom of each core photo- 
graph. Cores are arranged from proximal to the current source 
(core 5) to distal (cores 1 and 2; Fig. 1 and Supplementary 
material 1). Event Unit 1 is truncated in cores 3 and 5 based on 
depth of core penetration, and there true thickness is unknown 


% clay, % very fine silt, % fine silt, and % medium 
silt and the negative PC1 axis associated with above 
average charcoal, % coarse silt, % very fine sand, 
% medium sand, % coarse sand, and % very coarse 
sand. This dichotomy in sediment type suggests that 
two depositional environments dominate the sag pond 
record. The dominance of finer grain sizes associated 
with positive PC1 values suggests a low energy envi- 
ronment wherein sediment is deposited, accumulated, 
and preserved, likely in standing water or intermit- 
tent lake conditions. As a result, we interpret positive 
PC1 values as representative of ambient sedimen- 
tation associated with non-event processes such as 
playa or lacustrine deposition. For this paper, we are 
not concerned with the ambient sediments except for 
their value in focusing our selection of depths for the 
presence vs. absence of Erodium. Punctuating these 
periods of ambient sedimentation are several coarse- 
grained sediment units (negative PC1 values)—the 
focus of this paper. 

Coarser sediments dominate the negative PCl 
values suggesting a higher energy depositional envi- 
ronment such as shoreline migration via changes in 
water depth, precipitation-related runoff (i.e., event 
sedimentation), concentration of coarse sediment via 
deflation, or earthquake disturbances. Although the 
PCA identifies 7 units where the PC values are less 
than zero (Fig. 3), two of the negative PC1 units are 
markedly different—EU 1 and EU 2 (Figs. 3 and 5). 
EU 1 and EU 2 are noteworthy for several reasons: 
(1) they are the most negative PC1 intervals, (2) they 
are thick enough to extract meaningful sedimento- 
logical information (i.e.,> 1 cm or>one point of data 
per event unit), (3) they have erosional lower bounda- 
ries, (4) they contain the coarsest sediment, (5) they 
are normally graded, (6) they contain low % TOM 
and lower MS values, (7) they reveal an increase in 
charcoal near their upper boundary, and (8) they are 
visually correlated across all five cores at the study 
site. This combination of criteria suggests high 
energy deposition across the basin with diminish- 
ing energy from source (core 5) to sink (cores | and 


2; Fig. 1; Sabatier et al. 2022). We contend that the 
most likely explanation for EU 1 and EU 2 are precip- 
itation-related runoff events generated by SPEs. High 
intensity flow into the sag pond explains the erosional 
bases that characterize EU 1 and 2, as described for 
a mid-Holocene flood event in Lake Elsinore (Kirby 
et al. 2021). The thinning of the events from source 
(core 5) to sink (cores | and 2) indicates a progressive 
loss of energy as the runoff event flows down the long 
axis of the sag pond—a key criterium for discerning 
event origin (Sabatier et al. 2022). Normal grading 
is often associated with rapid deposition of sediment 
in aqueous environments wherein water—sediment 
energy dynamics generate grain-size separation (Wil- 
helm et al. 2012; Smith et al. 2019; Sabatier et al. 
2022). Evidence for laminae, especially in EU 2, sug- 
gest progressive energy loss and grain-size separation 
during sediment deposition (Sabatier et al. 2022). The 
presence of a thin clay layer in EU 1 suggests pulsed 
sedimentation, perhaps concurrent events, and favor 
subaqueous deposition allowing the finest sediment 
(i.e., clay) to settle and accumulate. The charcoal 
peaks near the top of EU 1 and 2 also suggest that 
the least dense materials (i.e., organic debris) were 
deposited as the flow ceased and the slow settling 
of low-density detritus commenced. The latter also 
favors deposition in a subaqueous environment. The 
large age range of charcoal dated in EU 2 suggests 
significant erosion of the drainage basin, tapping into 
older stratum containing older charcoal and/or ero- 
sion of old charcoal sitting on the landscape. Low 
organic contributions suggest dilution of organics via 
higher clastic sediment contributions during event 
deposition and/or post-depositional oxidation of the 
organics. Low MS values within the event units sug- 
gest a change in the mineralogy of event vs. ambient 
sedimentation and/or post-depositional alteration of 
magnetic minerals. Based on the above evidence, EU 
1 and EU 2 are interpreted as water-rich, sediment 
laden flows that likely generated turbidite-like units 
upon deposition, generated by precipitation-related 
runoff events, most likely into standing water (Fig. 5). 
However, we cannot unequivocally rule out that these 
event units were deposited on a playa surface as a 
surface flow event without the presence of significant 
standing water. In either scenario, the genesis remains 
linked to precipitation-related runoff in response to a 
SPE. 
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Fig. 3 Core CLPC21-4 with analyses versus depth. Key ages 
relevant to our discussion are shown to the far left in cali- 
brated years CE with the 95.4% confidence range. Presence 
or absence of Erodium is shown by green arrows (present) 
and red arrows (absent). °’Cs age is noted by black arrow; 
however, the Cs age is within an event unit is considered 
unreliable. All event units as identified by negative PC1 val- 
ues are highlighted in beige except for EU 1 (blue) and EU 2 
(green) =Event Unit 1 and Event Unit 2. From left to right: 


We also present, and refute, the three most rea- 
sonable alternative explanations for EU 1 and 2. The 
first alternative explanation for our two EUs is that 
they reflect a sedimentological response to sustained 
drought and the subsequent progradation of the shore 
zone across the core site (i.e., lake-depth shallowing; 
Kirby et al. 2018, 2023). If true, we expect to see a 
gradual change from fine grained clays and silts asso- 
ciated with a deeper water, low energy environment to 
coarser sandy silt and/or sand associated with a shal- 
lower water, high energy environment. As noted, the 
basal contacts for EU 1 and 2 are abrupt and erosive, 
and the units are normally graded (both units) with 
evidence for well-defined-to-faint laminae (only EU 
2). We contend that this abrupt and erosive transition 
to coarse sediment is not consistent with a sustained 
progradation of the shore zone basinward in response 
to a slow lake-level regression. As a small sag pond, 
changes in lake level are likely rapid—based on 
modern observations—and this rate of change will 
exceed the response/equilibrium time for establishing 
a wholly new grain size-energy depositional environ- 
ment. Moreover, the sag pond is extremely narrow 
across its long axis (<70 m) and likely would not 
record rapid and subtle changes in the position of 
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magnetic susceptibility (x 10’ m? kg.~!), percent total organic 
matter, charcoal counts per | g dry weight, percent clay, per- 
cent very fine (VFsilt) silt, percent fine (Fsilt) silt, percent 
medium (Msilt) silt, percent coarse (Csilt) silt, percent very 
fine (VFsand) sand, percent fine (Fsand) sand, percent medium 
(Msand) sand, percent coarse (Csand) sand, percent very 
coarse (VCsand) sand, and PC1 values. Grain size intervals 
according to the Wentworth scale (Wentworth 1922) 


the littoral zone, as the distance between the lake’s 
depocenter and shore is horizontally insignificant. 
On-site surface analysis as well as lidar imagery from 
the study site by Bevis and Hudnut (2005) reveal no 
observable beaches or stranded beaches to indicate 
prior, sustained highstands. Consequently, we do not 
think that EU | and 2 are shore-zone deposits associ- 
ated with a prolonged drought or lake-level lowstand. 
As a result of these observations and reasonings, we 
rule out shore-zone progradation as an explanation 
for the EU 1 and 2. 

Our second alternative explanation is that the 
event layers represent an encroachment of the fan- 
delta observed at the northwest side of the sag pond. 
Fan-deltas represent the distal portion of an alluvial 
fan where the fan meets and interacts with a standing 
body of water. They generally contain poorly sorted, 
coarse sediment, lack grading, and show both lateral 
and vertical heterogeneity (McPherson et al. 1987). 
We do not favor this origin for three reasons. First, 
EU 1 and EU 2 are well-sorted with normal grading, 
and they express lateral and vertical coherence across 
the basin. Second, fan-delta migration likely requires 
a significant period—decades to centuries—to form 
as well as a period of relatively persistent runoff to 
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move, establish, and deposit the fan-delta. Third, fan- 
deltas, by definition, require a body of standing water. 
During a period of sustained wetness (i.e., pluvial), as 
lake level rises, the fan-delta system will not prograde 
but retreat upland in response to the change in base- 
level (i.e., the lake-level transgression). For these rea- 
sons, we do not favor the interpretation that our event 
units (EU 1 and EU 2) represent fan-delta migration. 

Our third alternative explanation for EU 1 and 2 is 
a seismogenic origin, such as a sand blow (i.e., ver- 
tical sand migration via seismogenic liquefaction) 
or seismite (i.e., soft sediment-deformation features 
associated with seismogenic, subaqueous shaking; 
Sabatier et al. 2022). We contend that a seismic ori- 
gin is unlikely for several reasons. First, both EU 1 
and 2 lack typical soft sediment-deformation features 
traditionally associated with lacustrine seismites 
(Bowman et al. 2004; Montenat et al. 2007; Sabatier 
et al. 2022). On the contrary, they are both normally 
graded and retain evidence, such as laminae in EU 2, 
for grain-size separation associated with water-energy 
dynamics during the flow of a water-rich, sediment- 
laden medium. Second, both units thin from source to 
sink suggesting a progressive loss of energy that is not 
readily explained by seismogenic processes such as a 
seismite (Migowski et al. 2004; Sabatier et al. 2022). 
Third, a seismite caused by shaking of the basin will 
not necessarily introduce coarse grains as it will only 
rework that sediment already in the basin. Not only 
are EU 1 and 2 coarser grained than the ambient sedi- 
ment, but each unit becomes progressively finer and 
thinner from source (core 5) to sink (core 1). Fourth, 
EU | and 2 are characterized by abrupt, erosive lower 
boundaries, which cannot be easily explained via 
seismogenic processes, such as a sand blow. There- 
fore, we suggest that a seismic origin is very unlikely 
for EU | and EU 2. 


EU 1| and EU 2 timing and the Little Ice Age 


Initial depth selection for ‘tC age control preceded 
completion of the sediment analyses and was based 
on visual inspection of the core and initial charcoal 
counts. As a result, six of the seventeen radiocar- 
bon ages were extracted from ambient sediment 
units as defined by positive PC1 values (Table 2). 
To determine the age of EU 1, we focus on the ages 
collected from ambient sediment below the EU 1’s 
erosional base (Figs. 5 and 6). Charcoal from 154 


to 153 cm and combined 155 to 153 cm provide a 
range from 1470 to 1640 CE and 1510 to 1590 CE, 
respectively, at the 95.4% probability distribution 
under the calibration curve (Table 2). If the event’s 
erosional base reflects the removal of sediment 
during sediment emplacement, we consider these 
ages to reflect the event’s maximum age and likely 
younger when sediment erosion is considered. As a 
result, we assign a maximum age for EU 1 at 1470 
CE and a minimum age at 1640 CE. Based on this 
range, it is possible that EU | represents the terres- 
trial equivalent to the 1525 CE SBB flood event. 
Eleven of the seventeen C ages were subse- 
quently identified as occurring in event sediment 
units as defined by negative PC1 values (Table 2 and 
Fig. 6). The ages within event sediments are likely 
reworked and represent a variety of materials and 
their respective ages on the landscape at the time of 
erosion. Despite this reworking, we can infer impor- 
tant information about the event’s age by examin- 
ing the distribution of ages in the event unit. The 
youngest material dated within an event provides a 
maximum age for the event unit, assuming that it is 
more difficult to rework younger material into older 
sediments. Using this rationale, the youngest age in 
an event unit necessarily requires that the surface 
age of the landscape undergoing erosion was at least 
as old as the youngest date. For example, EU 2 con- 
tains five '4C ages from detrital charcoal or charred 
wood (Table 2 and Fig. 6). These ages (n=5) range 
from 1440 to 1800 CE (Table 2). As a result, we 
infer a maximum limiting age for EU 2 using the 
distribution of ages from the youngest sample. The 
youngest sample occurs at 129-127 cm with a mini- 
mum to maximum age range of 1740-1800 CE at 
the 95.4% probability distribution under the calibra- 
tion curve. Consequently, we suggest that EU 2 is 
likely no older than 1740 CE and no younger than 
1800 CE (Figs. 5 and 6). Moreover, our suggested 
age for EU 2 is independently supported by the 
first occurrence of Erodium at 108-107 cm ranging 
from 1750 to 1765 CE (Mensing and Byrne 1998). 
Assuming the Erodium age—as determined from 
SBB marine sediments—is reasonably applicable to 
the Carrizo Plains (100 km NNE of SBB), we can 
additionally constrain the age of EU 2 to 1740 to 
1765 CE. Whether the former (1740-1800 CE) or 
latter (1740-1765 CE), the inferred age of EU 2 
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<Fig. 4 Scatter plots comparing various sand sizes (including 
percent total sand) versus percent clay based on PC1 values. 
Negative PC1 values =event sedimentation (red circles); posi- 
tive PC1 values = ambient sedimentation (black circles). a per- 
cent clay versus percent very fine sand, b percent clay versus 
percent fine sand, ¢ percent clay versus percent medium sand, 
d percent clay versus percent coarse sand, e percent clay ver- 
sus percent very coarse sand, and f percent clay versus percent 
total sand 


is remarkably similar in age to the 1760 CE flood 
layer identified in the SBB (Sarno et al. 2020). 

Notably, EU 1 and EU 2 are separated by only 
2 cm of ambient sediment in core 4. Because EU 2 
is characterized by a sharp, erosional base, it can be 
inferred that some amount of sediment (i.e., time) 
was removed during event deposition and the 2 cm 
of remnant ambient sediment was originally thicker 
(i.e., more time between EU 1 and 2). Without know- 
ing the sedimentation rate and acknowledging the fact 
that our study cores are punctuated by event sedimen- 
tation, it is difficult to estimate the amount of time 
separating EU 1 and EU 2 and the amount of sedi- 
ment removed via erosion during EU 2. However, the 
thickness of the ambient sediment unit between EU 1 
and 2 increases as the core distance from the source 
lengthens, suggesting a reduction in erosive power 
from source to sink (core 5=proximal, core 1 and 
2=distal; Fig. 2). This observation provides addi- 
tional and independent support that EU 1 and EU 2 
are separated by a measurable period of time, and 
they are not genetically linked by a temporally cor- 
relative series of events. As a result, our best esti- 
mate for EU 1 is between 1470-1640 CE and EU 2 
between 1740 to 1765 CE. 

The Little Ice Age (LIA) was a period of climate 
change during the late Holocene between 1450-1850 
CE (IPCC age range, Masson-Delmotte et al. 2013), 
characterized by changes in temperatures, hydrocli- 
mates, and atmospheric/oceanic circulation (Haug 
et al. 2001; Chiang and Friedman 2012; Koutavas 
and Joanides 2012; Rustic et al. 2015; Cvijanovic 
et al. 2017; Loisel et al. 2017; Slawinska and Rob- 
ock 2018; Lapointe and Bradley 2021). In Southern 
California, several paleoclimatic archives character- 
ize the LIA as reflecting cooler and/or wetter hydro- 
climates. For example, lake sediments from Abbott 
Lake (1450-1570 CE, 1710-1775 CE; Hiner et al. 
2016), Zaca Lake (1450-1500 CE, 1580-1710 CE; 
Kirby et al. 2014), and Lake Elsinore (1450-1530 


CE, 1580-1610 CE, 1670-1700 CE, 1750-1780 CE; 
Kirby et al. 2010, 2018) suggest enhanced precipita- 
tion-related runoff during the LIA. At the same time, 
ephemeral lakes in the Mojave Desert re-emerged 
(e.g., Silver Lake and Cronese Lakes) and wetlands 
expanded (i.e., Soda Lake; ~ 1650 CE; Enzel et al. 
1992; Miller et al. 2010; Honke et al. 2019). East of 
the Mojave in Arizona, Ely (1997) observed a sharp 
increase in flood frequency after 1450 CE. Tree ring 
reconstructed water year precipitation for the San 
Gabriel drainage basin shows more frequent, above- 
average precipitation during the LIA (Meko et al. 
2017). Similarly, tree ring Palmer Drought Sever- 
ity Index (PDSI) reconstructions suggest sustained 
periods of less drought and a generally above-aver- 
age PDSI for Southern California, except between 
1580-1720 CE (MacDonald et al. 2008). In the SBB, 
Sarno et al. (2020) identified only two floods within 
the LIA (1450-1850 CE sensu stricto) at 1525 and 
1760 CE. The 1525 CE SBB flood-event layer is the 
second thickest event layer (out of 12) in the past 
2000 years (Sarno et al. 2020). 

Our research from the Carrizo Plains reveals two 
event units likely related to SPEs that are contem- 
poraneous with the LIA chronozone (i.e., IPCC age 
range, Masson-Delmotte et al. 2013). Without older 
sediment studies at our site, we cannot say with cer- 
tainty that the Carrizo Plains were wetter-than-aver- 
age during the LIA. However, based on the spatiotem- 
poral signature of a wetter-than-average LIA across 
Southern California, it is not unreasonable to assume 
that the Carrizo Plains were also wetter-than-average 
during the LIA. Future work on longer Carrizo Plains 
sediment cores is required to test the latter statement. 
As we note above, our data suggest another five PC1 
negative event units in the Carrizo Plains study site. 
Without additional age control per core and addi- 
tional in-depth analyses per core, we cannot—at this 
time—assign a definitive origin to these additional, 
younger five PC1 negative units. 


Broader implications 


The results of this study demonstrate the utility of 
sag ponds as possible paleoclimatic archives. These 
ponds may prove particularly useful for document- 
ing SPEs in the geologic record due to their focused 
drainage basins and their occurrences in regions 
devoid of more traditional lake environments. 
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Fig. 5 A close-up view of EU 1 and EU 2 in core CLPC21-4 
(split halves of the same core) with analyses versus depth and 
digital photograph (colors enhanced to highlight EU 1 and 2). 
Key ages relevant to our discussion are shown in calibrated 
years CE with the 95.4% confidence range. The best-estimate 
age range is shown for both EU 1 and 2 in the colored boxes. 


Coastal and marine archives, such as SBB, more 
likely integrate regional runoff events related to 
SPEs from throughout the Southern California 
region and its more extensive drainage basin. How- 
ever, given the narrow width of SPEs (specifically 
ARs), the SBB flood events do not necessarily rep- 
resent the whole region and highlight the necessity 
to locate and study the spatiotemporal signature 
of SPE floods recorded in terrestrial archives. Sag 
ponds have the potential to improve our understand- 
ing of past SPE storm tracks in addition to SPE 
storm frequency. Adding records from other sag 
ponds in this region will allow us to disentangle 
the spatial signature of SPEs across coastal Califor- 
nia. Sag ponds may prove a valuable and generally 
untapped paleo archive that can improve our ability 
to distinguish the spatiotemporal signatures of SPEs 
in the geologic past. 
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From left to right: magnetic susceptibility (x 1077 m? kg~4, 
percent total organic matter, charcoal counts per 1 g dry 
weight, percent clay, percent very fine sand, percent fine sand, 
percent medium sand, percent coarse sand, percent very coarse 
sand, and PC1 values. Grain size intervals according to the 
Wentworth scale (Wentworth 1922) 


Conclusions 


Our study using sediments from a small sag pond 
along the SAF provides evidence for at least two SPE 
flood units over the past 500 years, within the chrono- 
zone of the Little Ice Age. Alternative explanations 
such as shore-zone progradation in response to sus- 
tained drought, fan-delta migration, or seismogenic 
origins are not supported by the visual and sedimen- 
tological data. Age control as determined from both 
event sediments and ambient sediments was used to 
constrain the events’ timing. Based on ambient sedi- 
ments below the erosive base of EU 1, we assign a 
temporal range between 1470 and 1640 CE, possi- 
bly correlative to the 1525 CE SBB flood unit. Using 
the distribution of 4C ages from within EU 2 and 
the first occurrence of Erodium, we assign a range 
between 1740 and 1765 CE, possibly correlative to 
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OxCal v.4.4.4 Bronk Ramsey (2021); r:5 Atmospheric Data from Reimer et al. (2020) 
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Fig.6 All 'C ages calibrated to years CE using OxCal 
(v.4.4.4) with the 95.4% confidence range (black o= median, 
black+ =mean). First appearance of Erodium is shown with 
a red+and its approximate age range (Mensing et al. 1998). 


the 1760 CE SBB flood unit. The timing of these 
events within the LIA agrees with regional archives 
that suggest a wetter, stormier LIA hydroclimate. 


Ambient versus event sedimentation as determined by PC1 
values are shown along the left y-axis. Proposed maximum age 
of EU 1 and EU 2 are shown bolded text at the 95.4% confi- 
dence range 


However, without older sediments pre-dating the LIA 
at our site, we cannot say unequivocally that the LIA 
was wetter-than-average in the Carrizo Plains. This 
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study demonstrates that sag ponds represent a viable 
archive for capturing and preserving SPE flood units, 
especially in regions where natural lakes are rare— 
such as Southern California. Future work should fol- 
low a similar protocol as that used here—core tran- 
sects along the sag pond’s long axis, *C ages from 
ambient and event sediment, and varied sedimento- 
logical analyses. 
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